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We connect three phenomena of wave packet dynamics: Talbot images, revivals of a particle in 
a box and fractional revivals. The physical origin of these effects is deeply rooted in phase factors 
which are quadratic in the quantum number. We show that the characteristic structures in the time 
evolution of these systems allow us to factorize large integers. 
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I. TALBOT EFFECT, REVIVALS AND FACTORIZATION 

Although so much has been explained in optical science by the aid of the undulatory hypothesis, 
yet when any well-marked phasnomena occur which present unexpected peculiarities, it may be of 
importance to describe them, for the sake of comparison with the theory. 

This quotation is the opening sentence of a paper jlj by H.F. Talbot entitled Facts relating to optical sciences, 
published in 1836. In this article Talbot reports his experiments on interference of light. We quote again from 
his seminal paper: 

About ten or twenty feet from the radiant point, I placed in the path of the ray an equidistant 
grating* made by Fraunhofer, with its lines vertical. I then viewed the light which had passed through 
this grating with a lens of considerable magnifying power. The appearance was very curious, being 
a regular alternation of numerous lines or bands of red and green colour, having their direction 
parallel to the lines of the grating. On removing the lens a little further from the grating, the bands 
gradually changed their colours, and became alternately blue and yellow. When the lens was a little 
more removed, the bands again became red and green. And this change continued to take place for an 
indefinite number of times, as the distance between the lens and grating increased. In all cases the 
bands exhibited two complementary colours. It was very curious to observe that though the grating 
was greatly out of the focus of the lens, yet the appearance of the bands was perfectly distinct and 
well defined. 

* A plate of glass covered with gold-leaf, on which several hundred parallel lines are cut, in order to transmit the light 
at equal intervals. 

In today's language Talbot has considered the diffraction of light from a grating in the near-field zone. He 
found that the intensity distribution immediately after the grating repeats itself periodically at multiples of 
a characteristic distance, later called Talbot length. Moreover, at fractions of the Talbot length the pattern 
reappears in a rescaled version. Indeed, at the fraction 1/r of the Talbot length we find r substructures of the 
original pattern. The theoretical explanation of this phenomenon was provided in 1881 by Lord Rayleigh ||. 

The Talbot effect occurs not only for light waves, but also for matter waves. It has been observed for atoms 
j|, and most recently for Cgo molecules Q. For a comprehensive review of the Talbot effect, we refer to ref. ||. 

A phenomenon closely related to the Talbot effect occurs in the time evolution of wave packets |(| , such as a 
Rydberg electron. An initially well-localized packet spreads over its orbit but regains its original shape after a 
characteristic time, which is much larger than the classical period. Moreover, at fractions of this revival time, the 
wave packet splits into multiple copies. This effect of fractional revivals appears most clearly in the well-known 
problem of the particle in a box [Q |] . Revivals and oscillations can also be used for a new type of interferometer 
for light in a planar multimode waveguide |[|. 

Fractional Talbot images as well as the fractional revivals originate from phase factors which are quadratic in 
the quantum number. They give rise to Gauss sums |Tc| ], which have interesting number-theoretical properties. 

This immediately points to the newly emerging field of quantum information processing [ pd| , which has received 
an enormous drive from the discovery of the Shor algorithm [|l2| [l3j to factorize large numbers — a paradigm 
of number theory. The existence of distinct phenomena at fractions of a characteristic time suggests that the 
Talbot effect or wave packet physics may have links to the problem of factorizing numbers. Indeed, a recent 



proposal 1 14 1 makes use of a TV-slit Young interferometer, where N is the number to be factorized. Such a device 
is described by the Talbot effect due to a finite grating. The interference structure along the transverse direction 
at a fixed position after the grating serves as an indicator for a factor. For another approach we refer to ref. pi 
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In the present notes we show that the quadratic phase factors inherent in the time evolution of many quantum 
systems provide a tool to find the prime factors of a large number. Our article is organized as follows: In 
section II we briefly summarize the essential ingredients of the Talbot effect in the language of atom optics. We 
devote section 111 to the calculation of the free propagator of an array of wave packets. In section IV we connect 



this Green's function with the one of the particle in a box. In both cases we arrive at quadratic phase factors. We 
then show in section [v| that similar phases appear in many discrete quantum systems. They manifest themselves 
in the autocorrelation function. We dedicate section VI to cast the relevant sum into a form which brings out 
most clearly the phenomenon of fractional revivals. This form allows us in section VII to test our factorization 
scheme. We suggest a different approach towards factoring in section VIII based on properties of Gauss sums. 
We conclude in section IX by presenting a brief outlook. 



II. MODEL OF TALBOT EFFECT 



In order to set the stage for the mathematical treatment of the Talbot effect presented in the next section, 
we now define the principle setup. We consider the diffraction of an atomic wave from a grating. In principle 
our treatment is also correct for electromagnetic waves in paraxial approximation. Indeed, in this limit the 
d'Alembert wave equation reduces to the Schrodinger equation. For the influence of the higher order corrections 
we refer to ref. ]lr| . 

The grating, aligned along the s-axis, could be a mechanical or an optical one. Since we are interested in the 
subsequent propagation, it suffices to assume that the grating creates a wave function 4>(x) with the period d of 
the grating. In the case of an infinite array of slits, whose width is much smaller than their separation, the wave 
function 



f(x-nd) (I) 



consists of an infinite number of independent initial wave packets <p(x) separated by d. 

We assume that initially the atomic wave is under normal incidence, that is the wave vector fc, aligned along 
the z-azis, is orthogonal to the grating. When we consider an atom with an energy much larger than the recoil 
energy of the grating we can treat its motion along the z-axis classically. In this case time translates into the 
z-coordinate of the atom pJq , that is z = vt where v is its velocity along the z-direction. 

We now consider the distribution of atoms along the x-axis for fixed propagation time t. For a defined velocity 
v of the incoming atoms this time corresponds to a fixed position z behind the grating. At the Talbot time 
T, that is at the distance zt = vT the initial wave function <p(x) repeats itself. At time T/2 the interference 
pattern is identical to the initial one but shifted by half a period. At fractions of the Talbot time the period 
of the initial wave packet is a fraction of the original one. In the following sections we derive these results and 
show that they are a consequence of an intricate interference of phases. 

III. MATHEMATICS OF THE TALBOT EFFECT 

In the present section we briefly review the mathematical treatment of the Talbot effect. Our key tool is the 
Poisson summation formula. This analysis provides the foundation for the next section, where we emphasize the 
close connection between the Talbot effect and the particle in a box. 

A. Free time evolution of a periodic structure 

We consider the free motion of the wave packet 

oo 

<P(x-ruC) (2) 

1— — OG 

consisting of an array of infinitely many identical partial waves (f(x) separated by a distance d. The propagator 

G ficc (x, t\y, t = 0) = Af(t) e M*)(*-i/) a (3) 



(X) = 

n— — OG 
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of the free particle of mass M with the scaling factor 



S ^ (4) 



and the normalization 



allows us to find the wave function 

oo 

${x,t)= J dyGf toc (x,t\y,t = 0)<l>(y) (6) 

— OG 

at a later time t. 

Indeed, when we substitute the initial wave function, eq. (^|), and the Green's function, eq. (^), into the 
propagation equation, eq. we arrive at 

oo oo 
/, oo „ oo 

#r,t) = I AytfeW'-tf E V(y-nd)= I dyAf £ e^-f-^V^) (7) 



n— — oo , , n— — oo 



where in the last step we have introduced the integration variable y = y — nd. 
Therefore, eq. (0) takes the form 



ip(x,t) = J dyG Talhot (x,t\y,t = 0)<p(y) (8) 

— oo 

where 

oo 

G Talhot (x,t\y 7 t = 0)=N(t) e la ^-y- nd ^ (9) 

n— — oo 

denotes the Talbot propagator. 



B. Quadratic phase factors 

In order to bring out the relation of Gxaibot to the propagator Gbox of the particle in a box, discussed in the 
next section, we rewrite eq. (^|) with the help of the Poisson summation formula 

oo 

oo oo „ 

E /»= E / dn/He- 2 ™ (10) 

n— — oo m— — oo 

— oo 

where /(n) is a continuous extension of f n such that f(n) = /„ at integer values n. 

This relation allows us to replace the summation by a sum of Fourier integrals. Hence, the Green's function 
GTaibot, eq. (§), takes the form 

oo °° 

y.(x — y — nd) 2 ^ — 2-Kimn 



G Ta xbot(x,t\ y> t = Q) = E / dneM 

Til— — OO 

— OO 

oo 



where we have introduced the new integration variable £ = x — y — nd and the wave number 

k,„ = (12) 
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Moreover, we have recalled the definition, eq. (||), of the normalization Af(t) . 
When we perform the Gauss integral 



exp 



4a 



and make use of the definition, eq. (m, of the parameter a, we can identify the phase 



(7wc m ) 2 t _ t _ 



4a 2M h 



— = Em — = 2irm — . 



T 



Here we have introduced the energies 



and the Talbot time 



E,, 



2M 



T = 



Md 2 

hit 



Hence, the propagator GTaibot takes the form 

G Ta ibot (a;, t\y, t = 0) = - ^ ex P 



-im—{x — y) 



. , t 

exp [ — 27TITO — 



In the time dependent phase factor we recognize the quadratic dependence on the summation index m. 



(13) 



(14) 



(15) 



(16) 



(17) 



C. Integers and half integers of Talbot time 

According to eq. (|l7j ) the propagator for t = s ■ T with an integer s is identical to the propagator at time 
t = 0. Consequently, at integer multiples of the Talbot time the wave packet regains its initial shape. 
Furthermore, at t = (s + -g) • T the relation 



exp ^— 2 r Kim 2 {s + —) \ = exp (— 7rim 2 ) = (— l) m = exp (~imir) 
casts the propagator into the form 

^ oo 

t (a; > t=(s + l/2)r| 1 /,t=0) = - ^ exp 



G 



Talbot \ 



m— — oo 



2tt 



(18) 



(19) 



We can identify the right-hand side of this equation with the propagator at time t = shifted by half the period 
d. Hence, the initial wave packet repeats itself at half integer multiples of the Talbot time but is displaced by 
d/2. 



IV. PARTICLE IN A BOX 



We now turn to the standard problem of a particle of mass M in a box of length L. We first briefly review 
the essential equations and then derive the corresponding propagator. We conclude by comparing this result to 
the Talbot propagator. 
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A. Propagation 



The time evolution of an initial wave packet ip(x, t = 0) = f(x) in a box reads 

oc 

ik{x,t) =J2^nU n (x)e- iE - t / H , (20) 

n=l 

where the expansion coefficients 

L 

ip n = J dyip(y)u n (y) (21) 
o 

are in terms of the energy eigenfunctions 

u n (x) = y^sin(fc„a;) (22) 



with wave numbers 



and energies 



k„ =n~f (23) 

Li 



En = IF ' (24) 

When we substitute the expansion coefficients ip n , eq. (pl|), into the expression, eq. (p0|), for the wave function, 
we find the propagation equation 

L 

i){x, t) = Jdy G hox (x, t\y, t = 0) <p(y), (25) 
o 

in which 

2 °° 

G hox (x, t\y, t = 0) = j-J2 MKx) sm(k n y)e~ lE " t / h (26) 
n=l 

is the Green's function of the box. 

The quadratic dispersion relation, eq. (|2~i|), together with the definition, eq. (p3|), of the wave number, yields 
phases 

Ej- = 27rn 2 | (27) 



ith characteristic time 



4A1X 2 , x 



Hence, the Green's function of a particle in a box reads 

2 °° / t\ 

Gboxix, t\y, t = 0) = — 22 sin(fc n a;) sm(k n y) exp f -27rm 2 - ) . (29) 

n=l 



In complete analogy to the Talbot effect, we find that for t = T the phases, eq. (|27|), are integer multiples of 2n. 
At this time the Green's function is identical to the one at time t = 0. Consequently, the initial wave function 
revives, that is ij){x,t = T) — ip(x,t = 0) = <p(x). The similarity of the phases in the Talbot effect and the box 
problem justifies the use of the same symbol T for the Talbot time and the revival time. 

We conclude this section by noting that the quadratic dispersion relation, eq. ( |24|) , is also the origin of quantum 
carpets & 18 observed in a box (Tq, 19 and in many other quantum systems |20[ Ell. 



() 



B. Relation to the Talbot propagator 



It is instructive to cast the sum of products of sine functions in the propagator, eq. (|29|), into a slightly different 
form. For this purpose we recall the relation 

1 



sin(fc„:r) sm(k n y) 



e ik n (x-y) , e -ik n (x-y) _ e ik n (x+y) _ e -ik n (x+y) 



(30) 



The terms with negative phase can be combined with the ones having positive phase by extending the summation 
in eq. to negative values, that is 



1 x 

G ho x(x,t\y,t = 0) = — [ 



e ik„{x-y) _ e -ik n (x+y) 



n— — oo 



exp ( —2'ni'nf 



The term n = does not contribute since feo = 0- 
When we introduce the Green's function 



<>,xp j -2mn 2 ^ ) . 



(31) 



(32) 



G T (x,t\y,t = 0) = — ^2 ex P -inj(x-y) 

n— — oo 

eq. ([n]) takes on the form 

G hox (x, t\y, t = 0) = G T (x, t\y, t = 0) - G T (x, t\-y,t = 0). (33) 

Hence, the propagator of the box is the difference of two Green's functions with starting points y and — y, both 
leading to x at time t. This difference ensures the boundary condition that the wave function vanishes at the 
walls. 

When we compare the Green's function Gt, eq. (|3^), to the Talbot propagator Gxaibot, eq. (|l7|) , we find that 
they are identical provided L = d/2. 



V. TIME EVOLUTION AND AUTOCORRELATION FUNCTION 



So far we have concentrated on the time evolution of two specific quantum systems: a wave periodic in space 
and a particle confined to a box. Both systems have shown quadratic phase factors and complete revivals. In the 
present section we generalize this treatment and show that quantum systems with a discrete energy spectrum 
can display the same phenomena. 



A. Definition 



We consider a discrete superposition 

oo 

|V(* = 0)> =X>«l n > ( 34 ) 

71 = 

of energy eigenstates \n) with energy eigenvalues E n . Possible quantum systems are the inter-nuclear motion of 
a diatomic molecule p2| , a Rydberg electron [^3| , the center-of-mass motion of an atom in a standing light wave 
p4j or a single mode of an electromagnetic field in a cavity j25| . 
Due to the time evolution the energy eigenstates accumulate phases E n t/H leading to the state 



\m) = Y,^ e ' lEnt/r »- ( 35 ) 

71=0 

\s{t)\ 2 ^\m = Q)m))\ 2 (36) 

oo 

S(t) = W n e- lE ~ t/h (37) 

71 = 

is a measure for the overlap between the initial state \ip(t = 0)) and the state \ip(t)) at time t. This overlap is a 
sum over all quantum numbers n of weights W n = \ipn\ 2 with phases E n t/h. 



The autocorrelation function 



with 
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B. Quadratic expansion of the energy spectrum 



When the occupation probability W n has a dominant maximum around a quantum number n and the energy 
spectrum is only slightly changing with n, we can expand E n into a Taylor series 



E n £S E n 



On 



[n-n) + - -^ rr 
2 on z 



(n - n) 2 , 



retaining at most terms quadratic in n — n. With the definitions 



dE n 



and 



On 



1 d 2 E n 



Tel 



2 dn 2 



= n 



2n 



(38) 



(39) 



(40) 



where T c \ and T denote the classical period and the revival time, respectively, we arrive at the approximation 

(41) 



E n t/h ft; E n t/n + — (n - n) + — (n - n) 2 



of the phases. 

We substitute this expression into the definition, eq. (|37|), of S(t), introduce the new summation index m 
n — n and arrive at 



S(t)ne- iE ^ h S(t), 



where 



-2ni 



m m 



(42) 



(43) 



with W m = \tp m +n\ 2 ' Here we have extended the lower bound — n of the summation to — oo since the dominant 
contributions arise for m rj 0. 

In complete analogy to the Green's functions, eqs. (|l7| ) and (|3l]), of the Talbot effect or the particle in a box, 
the autocorrelation function, eq. (|3rj), is determined by a sum where the summation index enters the phase in 
a quadratic way. However, in contrast to these two examples, we now also have a linear contribution providing 
two distinct time scales T c \ and T. 



VI. FRACTIONAL REVIVALS 

In ref. |2(| we have devised a method to rewrite the sum S{t), eq. (]43|), in an exact way as to bring out the 
features of S(t) typical for the different time domains. We now concentrate on times 

t = I Tci + At = 1-T + e q/r T cl + At (44) 

that are close to a fraction q/r of T and are close to a large integer multiple £ of T c \. The contribution e q / r T c \ 
is a correction term, since in general we have 

lT cX +-T. (45) 
r 

According to ref. |^6| we can cast S{t) into the form 

oo 

S{t= q -T + e q/r T d + At)= Wn r Ml ] m (46) 
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with the Gauss sums 



r-l 



and the shape functions 

2&>(At)= / d/xW ? (/x)exp<J-27ri 



(47) 



At m\ / At\ T d 2 



(48) 



Here W(n) denotes the continuous extension of W m . 

The importance of this representation stands out most clearly for the example of a Gaussian weight function 



W(p) = 



1 



2vrAn 2 



exp 



2 VAn 



(49) 



of width An. 

In this case we can perform the integral, eq. (Ha), and find the explicit expression 



Z«(At)=^(At)exp 



(At - fT cl 
2a 2 (At) 



exp 



. (At - fT cl 
' 2a? (At) 



for the shape function. Here we have introduced the complex amplitude 

1 



Af(At) = 



y/l - i4irAn 2 (e q/r T c i + At)/T 



and the widths 



a 2 r (At) = 



1 



47r 2 An 2 



4An 2 



2 / e q/rT c \ + At 



T 



Ti 



and 



^(At) 



16ir 3 An 2 (e q / r T c i + At)/T 



1 a 2 e 9/'- T ci + At 
-An 1 — 

7T T 



Ti 



(50) 



(51) 



(52) 



(53) 



of the real and the imaginary Gaussians. 

According to eq. (|46|) the sum S(t), determining the autocorrelation function, contains the product 

W^X^At) of the Gauss sum Wm' and the shape function Tm (At). The latter consists of the product 
of a complex- valued square root, a real and an imaginary Gaussian. The Gaussians only take on non- vanishing 
values in the neighborhood of At = mT c \/r. When the separation T c \/r between two neighboring Gaussians 
is larger than their width ay they do not overlap. In this case the sum over m, that is, over the individual 
Gaussians, separates into a sequence of Gaussians. 

Hence, the autocorrelation function in the neigborhood of a fractional revival, that is at a time t = + 
e q / r T c \ + At consists of a sequence of Gaussians separated by T c \/r provided r is odd or 2T c \/r for r even. This 

dependence on the odd-even-property of r is a consequence B of the Gauss sums eq. @. When the 

Gaussians do not overlap the mth term in the summation, eq. (j4q), represents the mth fractional revival. 

When neighboring non- vanishing terms 2^ (At) and 2^/ (At) overlap, interferences between these terms arise. 
Then the phases of the complex Gaussian and the square root start to play an important role. Consequently, 
the sum S(t) exhibits a more complicated pattern. 



VII. FACTORIZATION USING WAVEPACKETS 



In the preceding section we have cast the sum S(t) determining the autocorrelation function into a sequence 
of complex-valued Gaussians. We now show, that this form suggests a scheme to factorize numbers. 
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A. General principle 

For this purpose we return to the non-overlap criterion and recall from eq. (|52]), that the width er r of each 
Gaussian is different since it depends on time At. The minimal width 

(mln) = Jcl (54 ) 



occurs for the time Ai m i n defined by 



-q/r 1 cl 



T cl + At mln = 0. (55) 



If £ q j r — the time Ai m i n of minimal width vanishes. Since the Gaussian with index m = also has its 
maximum at At = 0, this Gaussian is of minimal width. Moreover, it has no overlap with neighboring ones 
provided o> < -T c i. This condition puts the constraint 

An>^- (56) 



2tt 



on the width An of the Gaussian weight function W, eq. 
According to eq. (Ej), a vanishing correction term 



£q/r 

corresponds to the condition 



(57) 



I ■ r = q ■ N (58) 

where N = T/T c \ is the number to be factorized. 

Hence, a well-localized Gaussian at time £ ■ T c \ indicates e q / r — and thus provides the factor i of N. 

We conclude this section by briefly discussing the case when t is not a factor of N and therefore e q / r is 
non-zero. In this case the structure at At = is not a Gaussian of minimal width. Indeed, the widths of this 
Gaussian and its neighbors are so large, that they overlap considerably and the phase factors in N and in the 
imaginary Gaussian, eq. (|50|), lead to complicated interference structures. No clear fractional revival occurs at 
At = 0. 



B. Simulation 



In order to test our predictions, we encode the number iV in the ratio T/T c \ of the two time scales. We then 
search the autocorrelation function for maxima at integer multiples of T c j. 
In fig. |l| we show the autocorrelation function 

\Sn(t)\ 2 ^ 

as a function of dimensionless time t = t/T c \ for the Gaussian weight function, eq. (J49j). Here we have chosen 
the width An = 250. The number we want to factorize is N = 1309 = 7 • 11 • 17. Since the behavior of Sn(t) 
in the vicinity of an integer is important, we show \Sn(t)\ 2 around various integers. We recognize dominant 
maxima at r = 7, r = 11 and r = 17, which are indeed the factors of N. 



w ( m ) 

exp 



-2iri m H I r 

TV 



(59) 



VIII. GAUSS SUMS AND FACTORIZATION 



The success of our factorization scheme relies on the interference of quadratic phase factors. Why not concen- 
trate on the bare essentials of the method and eliminate the weight factors and the linear phase term altogether? 
We therefore consider the sum 



N-l 



SAr(n) = exp ^— 2-Kim 2 — ^ (60) 



m=0 
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FIG. 1: Autocorrelation function |<Sjv(r = I + Ar)| 2 , defined in eq. (|5S|), as a function of dimensionless time At in the 
vicinity of various integers I = 2, 3, 5 . . .. The goal is to find the factors of N = 1309 = 7 • 11 • 17. The autocorrelation 
function has a maximum at an integer, that is at the origin of each horizontal axis, provided this integer is a factor of N. 
In the present case we clearly recognize 7, 11 and 17 as factors. 

for fixed TV as a function of n. 

This sum was also investigated in ref. |27| and gives rise to the so-called curlicues. The emphasis of rcf. [[27i was 
on the self-similarity of the emerging structures. However, in the present discussion we focus on the possibility 
of finding factors by considering the real and imaginary parts of this sum. 

In fig. U we display real and imaginary parts of sn=21- We recognize that the factors 3 and 7 appear as the 
periods in the imaginary part. 

IX. CONCLUSIONS AND OUTLOOK 

Interference of quadratic phase factors emerging in the Talbot effect, the particle in a box or in fractional 
revivals of wave packets has the potential to factorize large numbers. So far our technique relies solely on 
interference and does not make use of entanglement, a purely quantum mechanical degree of freedom. In 



11 





FIG. 2: Absolute value of real part (left) and imaginary part (right) of the Gauss sum sjv(n), eq. (^0|), for TV = 21 = 3 ■ 7. 
The factors and their integer multiples appear in the imaginary part of S2i(n). 



a next step we want to combine the effect of the quadratic phases with the advantages entanglement can 
offer. We therefore have to consider composite quantum systems, such as highly dimensional spin systems, and 
develop generalized measurements (POVMs). Besides establishing a novel link to the problem of factorizing 
large numbers, this approach will povide new insight into the connections between quantum physics and number 
theory. 
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